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The programs in this report are for use 
within the Department of the Navy, and they are 
presented without representation or warranty of 



any kind. 



A PROCEDURE FOR ESTIMATING AN OBJECT’S POSITION 



BASED ON TWO OR MORE BEARINGS WITH A 
PROGRAM FOR A TI-59 CALCULATOR 

I . Introduction 

A procedure for estimating an object's position with 
bearings taken on or from two or more stations is developed 
in Section IV of this report. In the development of the pro- 
cedure/the following things are assumed: The object and the 

stations are fixed on the surface of a flat earth and the 
position of each station is known. The error in the bearing 
taken on or from a station is a normal random variable with 
a known standard deviation e and a mean of zero (if bias 
exists, it is known and removed); and station bearing errors 
are independent. The user instructions for a TI-59 program 
to implement the procedure are given in Section II, and the 
program listing is given in Section III. 

As an example to illustrate a use of the program, suppose 
bearings are taken on an object from three stations (1, 2 and 3) 
as illustrated in Figure 1. Also, suppose that the assumptions 
stated above are satisfied and that an initial estimate of the 
object's position is made and that it is relatively near the 
object. This assumption is discussed in Section IV. 



1 



Station 3 




FIGURE 1. Geometry for the Example 



Let the measured bearings and bearing errors (standard 
deviations) be: 



0 1 = 35° 

0 2 = 351° 
e 3 = 131° 

And let the ranges and bearing 




of the initial estimate be: 



2 



10,000 meters. 



38° 



r 

r 

r 



1 

2 

3 



15,000 meters 
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12,000 meters. 



3 2 = 346° 

B 3 = 127° . 



Use of the position estimation program with this data 
gives a final position estimate (fix) determined by: 

x = -512 meters 
y = - 75 meters 

where x is its East-West distance and y is its North-South 
distance from the initial position estimate. The East-West, 
North-South xy-coordinate system with its origin at the initial 
estimate is shown in Figure 2. So the final position estimate 
is 512 meters to the West and 75 meters to the South of the 
initial position estimate. 



North 

y 





Initial Estimate 


• Final Estimate 





FIGURE 2. The Location of the Final Position Estimate with 
Respect to the Initial Position Estimate. 
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Minimum area elliptical confidence regions for an object's 
position can also be found by using the TI-59 program. The 
centers of the regions are at the fix, and their axes lie along 
the x' and y' axes of the coordinate system obtained by 
rotating the East-West, North-South xy-coordinate system with 
origin at the fix through an angle Y . The angle Y is defined 
so that it is positive for a rotation in the counterclockwise 
direction . 

With the data from the above example, the program gives 
Y = -31°; so, the x' axis is directed 31° South of East. For 
a confidence region with minimum area and a confidence level of 
.9000, the ellipse bounding the region has a semi -major axis 
of 2064 meters, and a semi-minor axis of 1453 meters. The area 
of the region is 9.43 square killometers or 2.75 square nautical 
miles. The region is shown in Figure 3. 



North 




FIGURE 3. A .9000 Confidence Region for an Object's Position. 
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In the example discussed above, the position of the 
initial estimate is an input to the program. If this is 
not desirable, the program can be used to determine a po- 
sition for the initial estimate. The position is the in- 
tersection of the two bearing lines corresponding to the 
first two bearings entered in the program. Both options 
are illustrated in Section II. 

Since, in general, the smaller the bearing errors, 
the more likely that the initial estimate will be relatively 
near the object; small bearing errors can be considered to 
be a condition on the use of the procedure. 
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II. User Instructions 



The TI-59 program to which the user instructions in this 
section apply can be used to calculate the quantities described 
in Section I. 

The program requires the following inputs: 

1. the observed bearing from or on an object for two or more 
stations ; 

2. station positions relative to a reference position; and 

3. the bearing error (standard deviation) for each observed 
bearing . 

Station positions can be specified in either of two ways. 

In the first way, Mode A, each station's position is specified 
in terms of its bearing a and its range p from a reference 
position. In the second way. Mode B, each station's position 
is specified in terms of its East-West distance x (plus for 
East) and its North-South distance (plus for North) from a 
reference position. The reference position can be any convenient 
location. For example, if it were at a station, then for 
that station a = 0 and p = 0 or x = 0 and y = 0. 

The program also requires an initial estimate of the object's 
position. The user has two options: 

1. Let the program provide an estimate, or 

2. provide one with the input data. 

For Option 1, the initial estimate is at the intersection of the 
bearing lines determined by the first two observed bearings entered 
into the program. For this reason, if this option is chosen, the 
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first and second groups of data entered should correspond to 

the two stations estimated to have the smallest uroducts r.e.. 

1 1 

Although in this option the reference position cannot be at 
the initial estimate, it can be at one of the stations. If 
only two stations are involved, the final estimate is at the 
intersections of the bearing lines. (If the second option of 
either mode is used with an initial estimate which is not at 
the intersection of the two bearing lines, the coordinates of 
the final estimate will differ from coordinates of the inter- 
section to the degree of the approximations involved in the 
estimation procedure.) 

Two ways of providing confidence (probability) region data 
are available. In the first way. Mode C, the confidence 
(probability) p is specified. In the second way. Mode D, the 
multiplier k is specified where ko^, and ka~, are the 
semi-axes of the bounding ellipse. 

The values of various quantities calculated by the program 
are either stored in registers or appear in the display. If a 
PC-100A printer is used, some of these values will be printed. 

The location of calculated values and the printing format is 
given after the user instructions. Those quantities which are 
not described in Section I are described below in the User In- 
structions or in Section IV. 

All angles required or calculated by the program are 
in decimal degrees. 
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Step 



Instructions 



Enter Press Display 



1. If the calculator has been in use and 
flags have been set or the memory 
repartitioned, turn the calculator off 
and then on . 

2. Read Side 1 and Side 2 of Card 1. 

3. For confidence region calculations. 
Read Side 3 of Card 2. 



MODE A: Station Locations Specified in 

Terms of Bearing and Range from 
a Reference Point. 

4a. If the initial position estimate will 
be determined by the program, go to 





Step 


7a. 


See the 


note on Page 10. 






5a. 


Enter 


the 


initial 


estimate's bearing. 


a* 


A' 


6a. 


Enter 


the 


initial 


estimate's range. 


P* 


R/S 


7a. 


Enter 


the 


measured bearing on the 


6 . 

i 


A 



object from a station or the reciprocal 
of the measured bearing on a station 
from the object. 



8a. 


Enter 


the 


station ' s bearing . 


a i 


R/S 


9a. 


Enter 


the 


station's range. 


pi 


R/S 


10a. 


Enter 


the 


bearing error. 


e . 

l 


R/S 



11a. Repeat Steps 7a, 8a, 9a and 10a for all 
stations. The number of repetitions i 
appears in the display after Step 10a. 
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Step 



Enter Press 



Display 



MODE B: Station Locations Specified in Terms 

of East-West Distance and North-South 
Distance from a Reference Point. 



4b. 


If the initial position estimate will 
be determined by the program, go to 
Step 7b. See the note on Page 10. 






5b. 


Enter the initial estimate's East-West 
distance . 


X* 


B' 


6b . 


Enter the initial estimate's North- 
South distance. 


y* 


R/S 



7b. Enter the measured bearing on the 
object from a station or the 
reciprocal of the measured bearing 
on a station from the object. 

8b. Enter the station's East-West R/S 

distance . 

9b. Enter the station's North-South R/S 

distance . 

10b. Enter the bearing error. e^ R/S 

lib. Repeat Steps 7b, 8b, 9b and 10b for all 
stations. The number of repetitions i 
appear in the display after Step 10b. 



0 . 

i 



BOTH MODES 

12. Calculate the East-West distance, R/S 

the North-South distance, the 
bearing and the range of the 
position estimate relative to the 
reference position. Also calculate 
the rotation angle y , and the 

standard deviations C/s, and C/s. . 

x' y ' 

(See Section IV.) 
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Step 



Instructions 



Enter 



Press 



Display 



13. 

14. 

15. 

16. 

17. 

NOTE: 



For confidence (probability) region 
calculations, go to Step 14 if the 
confidence (probability) for the 
region is specified. If k is 
specified where kg-. and kg Ai 

X y ' 

are the semi-axes of the bounding ellipse 
with the larger the major axis, go 
to Step 16. 

Enter p, the confidence level p C Area 

(probability) and calculate k, 

kcu,, kn and the area of the region, 

x ' y 1 

(The area units correspond to the 
distance units used.) 

For a different value of p, go to 
Step 14. 

Enter k and calculate the confidence k D Area 

level (probability) p, koA, , kg,,. and 

x 1 y 1 

the area of the region. (The area 
units correspond to the distance units 
used. ) 

For a different value of k, go to 
Step 16. 



If a data entry error occurs in either mode, 
press RST and then use the following pro- 
cedure: For Option 1, return to Step 7 

and repeat all data entries. For Option 2, 
return to Step 5 and repeat all data entries. 

Also, if a position estimate is to be de- 
termined for a new object position, follow 
this procedure. 
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NOTES : 



a) The program printing format is given below: 

For the initial data, i = 1, 2, ... , n with one space 

between groups: 



Mode A 



a 5 



P’ 
0 . 

2 

a ■ 



e . 

l 



initial estimate 
if provided 



Mode B 



0 . 

l 



x . 

l 



y i 



The format for the calculated position data is: 



x 

y 



a 

P 

Y 
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For the confidence (probability) region portion of the program 
the format is after pressing either C or D: 



P 

k 



ka~ , 

x ' 


semi-axis 


<>1 

ID 


semi -axis 


Area 





b) The following data 


is stored in the indicated registers 


Data 


Registers 


X* 


R38 


y* 


R39 


Y 


R2 9 


X 


R30 


y 


R31 


a 


R32 


P 


R3 3 


a 


R16 


a 

y 


R17 



P 


R14 


k 


R15 


ka~ . 

x ' 


Rl 8 


ka~ 

y 


R19 

12 



R19 



Four data tapes for a sample problem are given below. 



Distance units have not been specified, but they could be 
meters for example. Angles are in degrees. Option 1 (initial 
estimate not provided) for Mode A and Mode B is indicated by 
A and by B and Option 2 (initial estimate provided) is indicated 
by A ' and B ' . 

For each mode and each option, the input data are 
indicated. The data determine the relative locations of three 
stations as well as the observed bearing of an object from 
each station. 

For A and B, the reference location is at Station 1 
and the initial position estimate (determined by the program) 
is at the intersection of the bearing lines for Station 1 and 
Station 2 . 

The intersection has coordinates x* = 906.4853528 
and y* = 17296.77092 with respect to Station 1. 

For A' and B', both the initial estimate and the 
reference location are at the intersection of the bearing lines, 
so a* = 0 and p* = 0 and x* = 0 and y* = 0. 

The data for A, B, A' and B' are all equivalent, 
and each solution gives the same data for a confidence (prob- 
ability) region calculation. A tape with confidence (prob- 
ability) region results for both Mode C and Mode D which 
correspond to A, B, A' and B' is given with the first four 
data tapes. 
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Following these data tapes, there is a data tape which 
illustrates the effects of using only the bearings for the 
first two stations. The tape is for Mode A, Option 1. The 
values stored in Registers 38 and 39 (the x and y 
coordinates of the intersection of the bearing lines from 
Station 1 and Station 2 with reference, in this case, to 
Station 1) are also listed on the tape (as well as given 
above.) And, as can be seen, the initial estimate and final 
estimate correspond. 
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A 



A' 




U. a* 

0. p* 
01 

i 83. ai 
17320. 50808 p x 



33. 02 

2 1 3 • ct 2 
20000. p 2 
3. e 2 

303. 0 3 

i29. 5887755 a 3 
8717.797886 p 3 



~ 3 • :! 2 . 8 3745 6 7 x 
8 3 4 . 05 8 6 8 3 5 y 
201.7584019 « 

898.0393111 P 

- 7. 325392259 y 
787. 3663757 a 
1 2 3 3 . 0 8 0 776 a 
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< X < >i 



B 



B' 



C or D 



. 


e i 


u . 


x i 


LI • 


*1 


4. 


e i 


33. 


®2 


9 9 8 6 . 285 3 4 y 


x 2 


523. 3595624 


^2 


3 • 


e 2 


303. 




7624. 94649 


X. 


1 1741. 38795 


^3 


O t 


e 3 


573. 5878927 


X 


16462. 71223 


y 


1. 995471723 


a 


16472. 70157 


P 


7. 325392245 


Y 


(■ v . “ 1 K b 7 * 1 * b 


a. 


1233. 080777 


a* 



y' 



y. x* 

U« y* 



-90b« 4853528 x x 
-17296. 77092 yi 

4. e x 

O . 0 2 

- 1 U 0 92. 7 o U 7 x 2 
~ 1 6 r‘ 7 o . 41 1 2 b y 9 

-i ^ 



303. 9 3 

6 r' 1 y . 4 6 i 1 3 7 X 3 
-5555. 382969 y 

3 
e 3 

-332. 8974589 x 
-834. 058690b y 
201. 7584019 a 
898. 0393184 p 



0. 9 P 

2, 145966026 k 
16 39. 661493 ko«, 
4 6 46. 14945 4 ko~ , 

1 4 U 4 6 3 64. 97 Area 

. 3646647168 p 

2 . k 

1574.732751 k 0 
2466. 161553 ko|! 
1 2 2 0 0517. 6 Area 



— "7* O cr -7* Q •7* *7* *7' u 
! . ^ -J C* 7 -2' ? 

787. 3663755 
1233. 080777 



Y 



a,. 



x 

y 



i 

» 



16 



A 



•— 1 B 


6 i 




906. 4853528 


R38 


u. 


“l 








0. 


Pi 




17896. 77092 


R39 


4. 


e l 








■j 

•J •— 1 ■ 


0 2 








• — i “? •“ i 

w i -j • 


a 2 








1 0000. 


P 2 








"l 

o. 


e 2 








9 Hi-,. 4 85 3528 


X 








17296. 77092 


y 








3. 


a 








17320. 508U8 


P 








20. 35750198 


y 








818. 886822 


a .v 








3092. 663848 


€ 









y' 



0. 9 


P 


2. i 45966026 


k 


1757. 3032 99 


ka. 


6 6 3 6 . 751549 


kcr? 


36639720. 91 


„ Y 
Area 


u ill b 4 b b 4 I* 1 1 b O 


P 


m 


k 


1637. '•73644 


k°- 


6185. 327696 


ka? ( 


8 1 824 8 5 i - ' . 22 


Ar&a 
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To obtain the results given in Section I, use Mode A 
and take the reference position at the initial estimate 
(a* = 0, p* = 0). Then a;L = 218°, a 2 = 166° and 
a _ = 307°. The data tape for the calculation is given below 
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A' 



0. 


a* 


0. 


p* 


35. 


e i 


2iS. 


a i 


i 0000. 


pi 


4. 


e i 


35 i. 


02 


166. 


Q.2 


15000. 


P2 


(' ft 


e 2 


131. 


0 3 


307. 


a 3 


12000. 


*3 


c 

■J. 


e 3 


-511. 96 i 356 


X 


-7 er a *7 cr 
1 w • *T •_= ! • J •_,» 0 


y 


261. 617739 


a 


51 r. 4 y y y 6 y r 


P 


31. 23492683 


Y 


677. 2632305 


0 . , 


3 6* 1 . 6 y y y 6* y 2 


X* 

a -. 




y 


0. 9 


p 


d . 145 'r : b b U d 6 


k 


1453. 383883 


ka A 


2063. 751628 


ko*' 


9422966. 381 


y 1 

Area 


. 86* i 6647168 


P 


< C. m 


k 


1 3 5 4 . 526461 


k a „ 


1923. 377726 


ko?' 



y i y •-+ 6* y 4 . 6< u 5 



Area 
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III. Program Listing 



Before entering the program, press 2nd and then CP or turn 
the calculator off and then on. Next enter 5 in the display, press 
2nd and then Op 17. This repartitions the calculator's memory so 
that the complete program can be entered. Two cards are required 
to record the complete program. The first card is for the position 
estimation portion of the program and for the first part of the 
confidence region portion of the program. The second card is for 
the second part of the confidence region program. 

Before recording the program, enter 6 in the display, 
press 2nd and then Op 17. This returns the calculator's memory 
to the normal partition (479.59). Returning the calculator to 
the normal partition allows the two program cards to be 
read in the normal partition without forcing. When the program 
is used, it repartitions the calculator so that Bank 3 registers 
are program registers. 
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000 


91 


R S 


050 


oo i 


: 6 


LBL 


051 


002 


18 


C 5 


052 


003 


69 


OP 


053 


004 


20 


20 


054 


005 




ST* 


055 


006 


00 


00 


056 


007 


92 


RTH 


057 


008 


76 


LBL 


058 


009 


19 


D 5 


059 


0 1 0 


69 


OP 


060 


Oil 


21 


21 


061 


0 1 2 




RC* 


062 


0 1 3 


01 


01 


063 


014 


92 


RTH 


064 


015 


/ b 


LBL 


065 


0 1 6 


10 


E 8 


066 


0 1 7 


65 




067 


0 1 8 


89 


11' 


068 


0 1 9 


55 


- 


069 


0 2 0 


01 


1 


070 


021 


1 J H 


:~s 


07 i 


U22 


00 


o 


072 


023 


95 


= 


i J 7 3 


024 


92 


RTH 


074 


025 


7 6 


LBL 


075 


026 


12 


B 


0 7 6 


027 


86 


STF 


077 


028 


01 


01 


078 


029 


i 


LBL 


079 


030 


1 1 


fi 


080 


031 


‘J ‘ j 

•J u— 


X i T 


08 1 


032 


fc* 


IFF 


082 


033 


00 


00 


083 


034 


00 


00 


084 


035 


56 


56 


085 


036 


05 


5 


086 


037 


69 


DP 


087 


038 


17 


17 


088 


039 


-7 

O i 


IFF 


089 


040 


03 


03 


090 


041 


00 


00 


09 1 


042 


54 


54 


092 


043 


47 


CHS 


093 


044 


09 


9 


094 


045 


42 


STO 


095 


046 


00 


00 


096 


U4 =‘ 


42 


STO 


097 


048 


01 


01 


098 


049 


02 




099 



STO 


too 


48 


EXC 


09 


101 


19 


1 9 


OP 


102 


•J V 


IHV 


O 


103 


97 


DSZ 


STF 


104 


09 


09 


00 


1 05 


01 


01 


X 1 T 


1 U b 


12 


12 


PRT 


107 


01 


1 


C 5 


108 


69 


DP 


R/S 


109 


21 


21 


PRT 


1 1 0 


9 8 


H D V 


C 8 


1 1 1 


91 


R/S 


R/S 


112 


42 


STO 


PRT 


113 


18 


18 


e * 


1 14 


65 




R/S 


1 15 


43 


RCL 


PRT 


116 


14 


14 


E 8 


117 


42 


STO 


C 8 


118 


20 


20 


IFF 


119 


•j o 


S I H 


03 


120 


75 


- 


01 


121 


43 


RCL 


64 


122 


10 


10 


D* 


123 


'“I 


IHV 


IFF 


124 


44 


SUM 


01 


125 


20 


2u 


00 


126 


m Z l o 


S I H 


88 


127 


65 


x 


- 


128 


43 


RCL 


I‘ 8 


129 


19 


19 


sr 


130 


95 


= 


+ / - 


131 


CT ” 

•J J 


- 


S I H 


132 


43 


RCL 




133 


20 


20 


D 8 


134 


A 


S I H 


GTO 


135 


95 


= 


00 


136 


42 


STO 


QQ 
»• _• 


137 


.j S ~s 


O C« 


+ ' 


138 


43 


RCL 


D 8 


139 


14 


14 


x ; t 


140 


39 


CDS 


n 8 


141 


65 




I H V 


142 


43 


RCL 


P/R 


143 


18 


18 


CE 


144 


75 


- 


= 


145 


43 


RCL 


COS 


1 4 6 


10 


10 


X 


147 


a y 


COS 


X i T 


148 


65 


V* 


— 


149 


43 


RCL 



42 

09 

69 

86 

00 

99 

I 8 

91 

99 

18 

91 

99 

18 

91 

99 

10 

18 

87 

03 

01 

64 

19 

87 

01 

00 

88 

“7 CT 

r • j 

19 

95 

94 

65 

19 

6 1 

OO 

99 

85 

19 

32 

19 

—• 

O l’ 

24 

95 

39 

65 

32 

95 



21 



150 


19 


19 


200 


22 


INV 


250 


49 


PRD 


151 


Of: 




201 


• 1 = :* 


P/R 


251 


26 


26 


152 




J- 


202 


42 


STO 


252 


49 


PRD 


153 


43 


RCL 


203 


19 


19 


253 


2 ? 


2 7 


154 


20 


20 


204 


10 


E 1 


254 


*- 5 u- 


X i T 


155 


. i~; 
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IV. A Development for the Procedure 



In the development for the estimation procedure given here, 
all angles are in radians and the assumptions stated in Section I 
apply. 

Figure 4 shows three bearing lines from the ith of n stations. 
One is the observed bearing line of an object. One of length goes 




FIGURE 4. Problem Geometry. 



to the origin of an xy-coordinate 
system located at the object's un- 
known position. And one of length 
r^ goes to an initial estimate 
with known position but unknown 
coordinates (x,y). Note, esti- 
mates for -x and -y estimate 
the object's position. To find 
estimates -x and -y, consider 
the arc coordinates u^ = A.^(0^-<t>^) 
of the observed bearing line and 
v^ = *.^(3 ^-<|k) of the bearing line 
to the point (x,y) . They are de- 
fined by the three bearing lines 
and the circle of radius which 

goes through the object's position 
and which is centered on the station 
as shown in Figure 4 . 
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By defining vr = ^(0^-$^) (all angles in radians), 

u = v. + w.. Note, 6.-8- is known, but is not. 

111 'll 11 

However, can be expressed in terms of x and y, and, to first 

order, = x cos B^ - y sin B^; so, if tan (Bj_ - <f>^) ~ for 

i = l,2,...,n, that is, if (x,y) is relatively near the object's 
position, u^ = ^(0^-B^) +x cos B^ - y sin B^ for i = l,2,...,n. 

In this development, the initial estimate is assumed to be 
relatively close enough to the object's position so that the above 
approximation for u^ can be used and so that fi ^ can be taken 
equal to r^. With this assumption, u^ is a function of: 0^, 

the observed value of a random quantity; the known parameters r^ 
and B^; and the unknown parameters x and y. 

If a distribution for the 0^ can be specified, then a 
distribution for the U. can be determined and standard estimates 

l 

A A 

x and y for x and y can be considered. Maximum likelihood 

estimates are discussed in this section. Each 0^ is taken to 

2 

be a normal random variable with mean A. and variance e. . And 

l l 

the n random variables 0^, i = 1,2,... ,n (one for each station) 

are taken to be independent. 

The likelihood for a sample 0 , , 0 ,, , ... , 0 is then 

^ 12 n 



L<e i' e 2 v - ; (^=— ) [- 1 1 (e i-v 2 / e i] 

and the likelihood for a corresponding sample u^ , , • • • , is 
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L (u, ,u 0 , . . . ,u ) 
12 n 



n 

n 

i=l 



( /2? a . 

1 





where a. = r.e. (with e. in radians) since u. = r . ( 0 . — <±> . ) . 

ill l 1111 

By definition, the maximum likelihood estimates of x 

/s. /\ 

and y are the estimates x and y which make L (u, , u„ , . . . , u ) 
2 2 1 2 n 

a maximum. In this case, making L (u, ,u~ , . . . , u ) a maximum is 

12 n 

equivalent to making I^(u?/o?) a minimum. So, to find x 

A 

and y, solve the following two equations for x and y: 



8 (fcn L) 
3x 



0 and 



3 Un L) 
3y 



The solutions are x = x and y = y, and x and y are the 
maximum likelihood estimates. With w. = r.(0.-3.) and 

i ill 

the conditions assumed above these two equations are linear 
equations in x and y. And, 



n 



1 [w^ + x cos 3^ - y sin 3^] (cos 3^)/o^ = 0 



and 



n 



[w. + x cos 3- - y sin 3-] (sin 3 . ) /o . = 0 

l l 1 l l l 



And, in terms of the following quantities: 

2 2 2 
A = Z (cos B^)/ck , B = E(sin 3^ cos 3 ^)/ck , 

C = Z (sin 2 3 i )/c? , D = E(w i cos 3 ^ /°\ < 

E = Z (w. sin 3 • ) /a? , 
i i ' i 
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the equations are : 

A A 

Ax - By = -D 
Bx - Cy = -E . 

So the solutions are: 

x = (BE - CD)/ (AC - B 2 ) 
y = (AE - BD)/ (AC - B 2 ) . 

A confidence region can be constructed about an estimated 
position. In order to indicate how this is done, a probability 
region about the true position will be considered first. 

/X A 

Note, x and y are values of random variables. If a 
new set of bearings 0^,0 2 ,...,0 n is observed (for a fixed initial 
estimate and object), in general, a new pair of values x and 
will be obtained. 

A A 

If X and Y represent these random variables, then 

~ I n ? 

X = £ (VI. /o.) (B sin 3 . - C cos 3 -) 

(AC-B ) 1 1 1 1 1 



/v 1 9 

Y = 5— £ (VI. /ai) (A sin 3 . - B cos 3 -) 

(AC-B ) 1 1 1 1 1 

with = r^(0^-3^). (W^ is the random distance intercepted 

along the ith arc between the bearing lines defined by 0^ and 
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Note, X and Y have a bivariate normal distribution, 

since they are a linear combination of the n normal random 

variables w i ,W 2 ' ’ * * ' W n ' ° r ec 3 u; *- valentl y of the n normal 

random variables 0 ,0 . Also E(W.) = r . (<J> . - 3 . ) . 

i ^ n 1 ill 

If 3^ = <f>^ for i = l,2,...,n, that is, if the initial 
estimate of the object's position is at the object's position, 

E(W^) = 0 for i = l,2,...,n. And, therefore, E(X) = 0 and 
E(Y) =0. So, in this case, the "location" of the bivariate 

A A 

normal distribution of a point (x, Y) , the random 
coordinates of the object's estimated position, is the same as 
that for the point (-X, -Y) and both are centered on the 

object's position. However, the "location" of the distribution 

/N /N 

of (-X,-Y) is independent of the location of the initial 

/N A 

estimate when the coordinates (-X,-Y) refer to a coordinate 
system with origin at the initial estimate. This fact simplifies 
the establishment of a confidence region about the location of 
an estimated position. 

A region of minimum area for a given probability of con- 
tainment of an estimated position can be determined. The region 
is bounded by an ellipse which is centered on the object's position 
and whose axes lie along the axes of an x ' y ' -coordinate system 
obtained by rotating the xy-coordinate system centered on the 

object's position through an angle y . In this system, a ^ = 0, 

„ . xy 

that is, X' and Y' are independent normal random variables. 
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The two coordinate systems are illustrated in Figure 5 . 
The coordinates of a point in the two systems are related by 
y' y x' = x cos y + y sin y 

y' = -x sin y + y cos y 

x ' 

X 

These relations, along with 

o * /s =0, imply: 

x ' y ' 

FIGURE 5. Rotation Geometry. 




a 



2 



x' 



2 

a 

/N 

y' 



2 2 2 2 
o cos y + 2a cos y sin y + a sin y , 

x xy y 



2 

2 2 2 
sin y - 2a cos y sin y + a cos y 

A I A A ' 1 A 1 



X 



xy 



and 



tan 2y 



2a 



A /\ 




A 



X 



2£Z 

2 

- a 

A 

y 



where y, the angle of rotation of the coordinate axes, is 
positive in the counterclockwise direction. 

With the initial estimate of the object's position at 
the object's position (6^ = <Jk, i = l,2,...,n), so E(W^) = 0 
and Var(W^) = a?, 
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and 

O /N /• N 

xy 

Using 



~ 2 ' 
X (AC-B ) 



„2 _ A 

° * 2 ’ 

y (AC-B ) 



B 

xy (AC-B 2 ) 



tan 2y = 2B/(C-A) for 3^ = 4>^ / i = l,2,...,n. 

With the object's position known and, hence, 4 k known 

2 2 

i = l, 2 ,...,n, the above equations for a , o , a and y 

A /\ A A 

x y xy 

be used, since the initial estimate of the object's position 

be taken as the object's position. 

With values for a^, a^, a^ and Y, values for 

x y xy x ' 

can be found by using the equations in the middle of 

y' 

Page 30. And then, the probability that an estimated 

position will be within an ellipse of semiaxes ko^ and ka^ 

x' y' 



and 

So, 

for 

can 

can 

and 



a 2 = - 
x (AC 



1 r 2 2 

5 -^- l (l/o") (B sin 3 . - C cos 3 . ) ^ , 

-B^V 1 1 1 



a 2 = 



y (ac 



1 n 2 

— 5 — j 1 (1/a.) (A sin 3. - B cos 3.) 
-B*) 1 1 1 1 



1 n 2 

— 5 — £ (1/a.) (B sin 3.- C cos 3 • ) (A sin 3-- Bcos3.) 
-B ) 1 1 1 1 1 



(AC-B ) 



the definition for A, B and C, the above be 



come 
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which is centered on the object's position can be found. It is 
2 

1 - exp(-k /2) . (This result follows from integrating the bi- 
variate normal density over the ellipse.) And the area of the 

ellipse is 7 rk 2 cu a* . 

x ' y ' 

/\ A 

Given estimates x and y found by using the relations 

on Page 28, an ellipse with semi-axes ko„ and kg^ centered 

x' y' 

on the point with coordinates (-x,-y) in a coordinate system 

with origin at the initial estimate and oriented as indicated 

2 

by y is a 1 - exp(-k /2 ) confidence region. This follows 

/N /\ 

from the bivariate normal distribution of -X and -Y which 

in this system is centered on the object’s position. The ellipse 

2 2 

is defined if g^ and are known (the covariance matrix 

x y xy 

is known). And to the degree of the approximations involved, 

this can be assumed to be the case. In particular, by assuming 

the initial estimate of the object's position is at the object's 

position, which is consistent with assuming (6 ^-<Jk) is small, 

2 2 

values for , o ^ and Y can be obtained by using 

x y xy 

the relations on Page 31. These values can then be used to 
2 2 

determine and o A by using the relations on Page 30. 

x ' y ' 

And, then, with a value for k, a confidence region can be 
constructed. To the degree of the approximations involved, the 
shape of the confidence region is independent of both the object's 
position and of the initial estimate of the object's position. 
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For the case where bearings are taken from the object 



on two or more stations, 0 is the reciprocal of the 
bearing taken from the object. 

A discussion for this and for other bearings only 
position estimation procedures for situations similar to the 
one considered here is given in Reference 1 listed below. 
Reference 2 gives an equivalent bearings only procedure. It 
also gives a range only procedure, a range and bearing procedure 
and HP-9830A programs with which to implement the procedures. 
Using the fix determined by two lines of bearing as the initial 
estimate was suggested by this reference. 



The equations used in the program to determine (x*, y*) , 
the coordinates of the fix, are: 



x* sin ( 02 ~ 0 ^) 


= [P-L 


sin 


r— 1 

CD 

1 

i— 1 

a 


sin 


0 2 




" tp 2 


sin 


^ a 2 -0 2^ ^ 


sin 


6 1 


y* sin ( 6 2 — ) 


r— 1 

CL 

II 


sin 


(a 1 -e 1 ) ] 


cos 


0 2 




- I«2 


sin 


(a 2 -e 2 ) ] 


cos 


6 1 
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